Peak Effect and the Transition from Elastic to Plastic Depinning 
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We demonstrate for the first time an observation of the peak effect in simulations of magnetic 
vortices in a superconductor. The shear modulus C66 of the vortex lattice is tuned by adding a fictitious 
attractive short range potential to the usual long-range repulsion between vortices. The peak effect 
is found to be most pronounced in low densities of pinning centers, and is always associated with a 
transition from elastic to plastic depinning. The simulations suggest in some situations that over a 
range of values of cm the production of lattice defects by a driving force enhances the pinning of the 
lattice. 
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One of the most important aspects of superconducting 
materials is their behavior in magnetic fields. In Type 
II superconductors, magnetic fields induce the forma- 
tion of quantized vortices through which magnetic flux 
may penetrate the system. The application of a cur- 
rent in the presence of vortices generates an effective 
force that causes them to flow through the supercon- 
ductor, thereby dissipating energy and spoiling perfect 
conductivity. Much work, both theoretical and experi- 
mental, has been devoted over the years to finding and 
understanding mechanisms that pin the vortices, so that 
the superconductor will remain dissipationless even in 
the presence of the field |Q . In the presence of pinning 
centers, the superconductor displays a critical current j c 
above which dissipation sets in; physically this current is 
proportional to a critical force F p at which the vortices 
become depinned and can move through the system. 

The competition between intervortex interactions and 
pinning by disorder results in a surprising and long stud- 
ied phenomenon known as the "peak effect" . As the crit- 
ical field or critical temperature at which a sample loses 
its superconducting properties is approached, in many 
situations one observes an enhancement of the supercon- 
ductivity just before it is completely suppressed. Thus, 
j c exhibits a peak as a function of field or temperature, 
just before it vanishes. The earliest understanding of this 
phenomenon, the collective pinning theory, involves the 
softening of the elastic moduli of the vortex lattice as the 
superconducting order is suppressed ||||] , so that the vor- 
tices may settle more deeply into the pinning potential 
and thus become more difficult to dislodge. One diffi- 
culty with collective pinning is that tearing of the lattice 
is ignored in estimating j c : i.e., the lattice depins elasti- 
cally, not plastically. In recent years this assumption has 
been increasingly questioned in the peak effect regime, 
particularly for high-purity superconductors (e.g., NbSe2 
and systems with strong pinning centers Al- 
though there is accumulating experimental evidence of 
plastic motion in the peak effect regime, its precise effect 
on the size of j c is not known, and there is disagreement 



as to whether j c is enhanced or suppressed [|]JT^] 

by the onset of plastic motion. 

In principle much of this debate could be settled by 
direct imaging of the vortices near the depinning crit- 
ical current. However, in the peak effect regime such 
experiments are exceedingly difficult because the order 
parameter is suppressed near the critical temperature or 
critical magnetic field. Numerical simulations thus of- 
fer a unique window through which one may view the 
qualitative behavior of the vortices [jl2] [l6| . In this work, 
we present results demonstrating for the first time (al- 
beit, in a two-dimensional geometry) the peak effect in 
a simulated vortex system, show conclusively that it is 
associated with a crossover from elastic to plastic mo- 
tion, and find that under different circumstances lattice 
tearing may enhance or suppress j c . 

In recent years, the peak effect has been associated 
with the proximity of the vortex lattice to a melting tran- 
sition |I(|[l7],|l8| . Direct simulations of depinning in this 
situation pose enormous practical problems because near 
a critical point one inevitably has large thermal fluctua- 
tions. To circumvent this problem, we take note that in 
nearly every theoretical approach to the peak effect, it is 
not actually melting itself but rather the softening of cqq 
and other elastic moduli near the melting transition that 
is responsible for the effect. We thus consider a system 
of vortices in which the interaction may be varied so that 
the elastic properties of the system may be tuned directly, 
without the introduction of critical fluctuations. To sim- 
ulate a large number of vortices (iVy — 1600) we confine 
ourselves to two-dimensional systems, so that our simula- 
tions are most directly applicable to very thin films or su- 
perconductors consisting of effectively decoupled layers. 
In contrast to previous simulations, we focus on pinning 
centers that are dilute compared to the vortex density, 
which is most appropriate for systems with strong pin- 
ning centers Q§Qg. 

The precise form of the intervortex interaction we use is 
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H int = -t J^y\n\R-R'\+A v e-\ R - A '\ 2 ^} (1) 

R^R' 

where R are the position vectors of vortices, e 2 is the 
strength of the logarithmic interaction, and A v and 
are the strength and the range of a short range attrac- 
tive interaction. For large enough values of A v the in- 
teraction may in principle be attractive over a range of 
vortex separations; however, in all the simulations we re- 
port here, A v is small enough that the net interaction 
is repulsive at all distances. A uniform background is 
assumed to cancel out the diverging energy due to the 
logarithmic interaction. Because of the long-range po- 
tential, the bulk modulus is formally divergent, while the 
shear modulus cqq may be shown to have the form ]2C| ] 

C66 = ^ { ^-^E^(#) 4 -(|)V^} (2) 

R 

where n v = 2/(v / 3a 2 ,) is the density of vortices. We take 
£[, = 0.5ao so that the shear modulus has values egg = 
(n„/8)(e 2 — 1.768^), and tune cqq by varying A v . The 
vortices interact with the pinning centers through the 
potential: 

H pm = -A p ^e-^ 2 ^ (3) 

R,r 

where A p is the strength of the pinning centers and r are 
the positions of the pinning centers [^lj . 

Systems with Ny = 900 and Ny — 1600 vortices and 
different numbers of pinning centers located randomly 
are studied by a simulated annealing molecular dynam- 
ics (MD) method. Periodic boundary conditions are im- 
posed and an Ewald sum technique |^o| is used to com- 
pute the forces and energies. To equilibrate the system, 
the temperature is lowered from above the melting tem- 
perature to kT = O.OOle 2 in about 30 consecutive steps 
through typically 10 5 MD steps. The depinning force 
is then measured at very low temperature, using a qua- 
sistatic technique as follows. The center of mass 

of the system is shifted in steps of O.Olao by imposing a 
driving force. At each step, the driving force is allowed 
to fluctuate while the center of mass is fixed, and 1000 
MD steps are allowed to pass to equilibrate the shifted 
system. The average driving force required to hold the 
center of mass at this position is then measured over 200 
MD steps. The average driving force increases approxi- 
mately linearly with center of mass shift until the system 
finds a new minimum energy configuration, at which time 
the required driving force drops sharply. The depinning 
force is then defined as the peak driving force observed 
in this process. 

Fig. ^ illustrates our results for the depinning thresh- 
old force F p as a function of the shear modulus cqq for 



several different choices of £ p , A p , and N p , the number of 
pinning centers. In our simulations, we find that the peak 
effect is most pronounced when the density of pinning 
centers is small compared to the number of vortices, so 
we focus our attention on simulations with N p = 50 and 
100. For increasing values of cqq, we expect the number of 
pinned vortices to decrease, as in collective pinning. This 
general trend is confirmed in the inset of Fig. 1, which il- 
lustrates the fraction of occupied pinning sites P ccupied> 
defined as the fraction of sites for which a vortex may 
be found within £ p /y/2 of the pinning site center. When 
the lattice depins elastically (i.e., when tearing and de- 
fect formation may be ignored) , one expects the threshold 
depinning force to be proportional to the density of occu- 
pied pinning sites and the maximum pinning force that a 
site may exert, f max = (^4p/^p)v / 2e -1 / 2 . In the main part 
of Fig. [l], f max has been scaled out of F p , and two curves 
proportional to Poccupiod are plotted for the N p = 50, 100 
data with the proportionality constants chosen to match 
the threshold force for the largest values of cee, where the 
depinning is most elastic. As may be seen, the pinning 
force matches the expectations for elastic depinning rea- 
sonably well down to cq§/ 'A p n v « 2.0. For N p = 100, F p 
decreases monotonically to zero as cqq — > 0, whereas for 
A p = 50 there is a clear tendency for F p to overshoot the 
elastic depinning estimate before dropping to zero. This 
non-monotonic behavior of F p vs. cqq is the peak effect, 
and one of the surprising results of this study is that it is 
more pronounced for lower densities of pinning centers. 
This latter result is in qualitative agreement with experi- 
ment, for which samples that are more weakly disordered 
exhibit stronger peak effects ||. 

The deviations from the elastic depinning estimate oc- 
cur because for small values of cqq , the lattice easily de- 
forms and tears, allowing motion without forcing out all 
the vortices trapped in pinning sites. Fig. illustrates 
the trajectories of the vortices for large, intermediate, 
and small values of cqq. For the largest values, the vor- 
tex lattice largely retains its order as it depins. As the 
maximum of the peak in F p is approached, a crossover 
from elastic to more plastic motion is observed, in which 
dynamically changing channels form where vortex motion 
takes place. The width of these channels decrease with 
decreasing cqq. This behavior is reminiscent of plastic 
motion observed in superconductors with ordered arrays 
of pinning center s [pf and in simulations of disordered 
Wigner crystals |23| . As the falling edge of the peak 
effect is entered, a new qualitative behavior emerges in 
which the active channels of motion of the vortices are 
no longer dynamic, and the motion becomes very much 
like river flow 

Scenarios in which the peak effect is associated with 
a crossover from elastic to plastic motion have been ad- 
vocated by several groups in the last few years p|,p|,^0[ . 
The present simulations strongly support this viewpoint, 
although the precise evolution of the flows with cqq dif- 
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fers in some important aspects from what previously has 
been supposed. In particular the onset of plastic mo- 
tion in the peak effect regime has been thought to be 
associated with either a monotonically increasing |p[JlO| 
or decreasing |||ll| F p with decreasing cee- Our simu- 
lations demonstrate that in a sense both scenarios are 
true. Plastic flow, when it first sets in with decreas- 
ing C66, is associated with an increasing critical current. 
This is particularly true for very low densities of pin- 
ning centers, for which F p is enhanced by tearing [ p5[ . 
For low enough values of cqq, however, river flow motion 
sets in and F p becomes proportional to cqq. The latter 
behavior is quite sensible once one recognizes from the 
simulations that the motion of vortices for the smallest 
values of cqq correspond to river flow through channels 
that do not change dynamically. In this case the depin- 
ning force comes about due to interactions of the "rivers" 
(moving vortices) with the "river banks" (stationary vor- 
tices), whose ability to hold the rivers in place decreases 
with decreasing shear modulus. 

A useful way of characterizing the depinning force for 
the smallest values of cqq (the "static river" limit) is to 
assume that if few of the vortices trapped in pinning sites 
are pulled free by the depinning driving force, then the 
only relevant length scale in the system at the depin- 
ning transition is d cx (n^/n^) 1 / 2 , the average distance 
between pinning sites. In particular this implies a char- 
acteristic displacement scale for the lattice in the pres- 
ence of the driving force u max = r c d, with r c a unit- 
less constant, above which defects are produced so that 
the lattice becomes depinned. The work done by the 
driving force must provide the elastic deformation en- 
ergy just before defects and depinning set in, so that 
F p n v r c d ~ c 66 r 2 d 2 , or F p n v cx c 66 (n v /n p ) 1/2 . Fig. 3 il- 
lustrates this scaling relation, and one may clearly see 
a collapse of the data onto a single straight line for 
C6e/A p n v < 1.0. The collapse of the data indicate that d 
is indeed the only relevant length scale in the static river 
flow regime. We note finally that in our limit of dilute 
pinning centers, river flow motion is possible for an un- 
melted vortex system (cqq > 0), in contrast to what has 
been speculated for systems with dense pinning centers 
§1- 

In summary we have reported the first simulations of 
the peak effect in a vortex lattice. We observe a peak 
in the depinning force near the smallest values of Cqq 1 
demonstrate with particle trajectories that this peak is 
associated with a crossover from elastic to plastic motion, 
and find that the peak is most pronounced for low pinning 
center densities. 
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03814 and a Cottrell Scholar Award of Research Cor- 
poration (HAF), and by Korea Science and Engineering 
Foundation 961-0202-008-2 and BSRI-97-2448 (MCC). 

Note added: After the submission of this work, a pub- 



lication |26| appeared reporting experiments on Nb us- 
ing neutron scattering to measure correlation lengths of a 
vortex lattice. It was found that in the peak effect regime, 
the correlation length corresponding to shear displace- 
ments decreases monotonically to a minimum through 
the rising edge of the peak effect. This observation cor- 
roborates our finding that the softening of cqq may be 
the controlling parameter in the peak effect. 
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FIG. 1. Depinning force F p as a function of shear modu- 
lus cm for Nv vortices and N p pinning centers of strength A p 
and range £ p . The maximum possible force f max and the num- 
ber of pinning sites per vortex N p /Nv are scaled out so that 
data should collapse onto a single curve if the lattice depins 
elastically. Solid and dotted lines illustrate the expected be- 
havior for elastic depinning. Inset: Fraction of occupied pin- 
ning sites Poccupicd for the groundstate configurations found 
by simulated annealing. 



FIG. 2. Trajectory plots for depinned vortices at dif- 
ferent values of C66, illustrating the evolution from elastic 
to plastic motion. Crosses represent the locations of pin- 
ning centers. A p = 0.03e 2 , N v = 1600, and iV p = 50, 
and (a) cm/{A p n v ) = 4.17; (b) c m /(A p n v ) = 1.96; (c) 
c 6 a/(A p n v ) = 0.48. 



0.4 



z 



• • • • □ 
■ ■ • 



0.0 1.0 2.0 3.0 4.0 

c 66 /(A p n v ) 



FIG. 3. Depinning force F p times inter-pinning-center 
spacing as a function of shear modulus cm- The symbols rep- 
resent parameters as used in Fig. 1. The figure clearly shows 
that in the region where cm is small F p is proportional to 
cm, and that the average distance between pinning centers is 
the only relevant length scale for "static river flow" depinning 
(see text) . 
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